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Abstract 

We prove Nishida's 1971 conjecture stating that almost all low-energetic motions of 
the anharmonic Fermi-Pasta-Ulam lattice with fixed endpoints are quasi-periodic. 
The proof is based on the formal computations of Nishida, the KAM theorem, 
discrete symmetry considerations and an algebraic trick that considerably simplifies 
earlier results. 



1 Introduction 

The Fermi-Pasta-Ulam (FPU) lattice is the famous discrete model for a continuous 
nonlinear string, introduced by E. Fermi, J. Pasta and S. Ulam (7j. It consists of a 
number of equal point masses that nonlinearly interact with their nearest neighbors. 
The physical variables of the FPU lattice are the positions qj of the particles, see 
Figure and their conjugate momenta pj. 




Figure 1: Schematic picture of the FPU lattice. 

Fermi, Pasta and Ulam were interested in the statistical properties of the nonlinear 
FPU lattice. In fact, they expected that it would attain a thermal equilibrium, as 
was predicted by laws in statistical mechanics. This means that the initial energy of 
the lattice would be redistributed and, averaged over time, equipartitioned among 
all the Fourier modes of the lattice, see [14] . They performed a numerical experi- 
ment to investigate how and at what time-scale this would occur. The astonishing 
result of their integrations was that there was no sign of energy equipartition at all, 
see (7j and [Tl|: energy that was initially put in one Fourier mode, was shared by 
only a few other modes. Moreover, within a rather short time all the energy in the 
system returned to the initial mode. This recurrent behaviour has been observed 
in experiments on the FPU lattice with quite small as well as very large numbers 
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of particles, on short and long time-scales, and we are led to believe that at low 
energy the FPU lattice behaves more or less quasi-periodically. This observation 
was a big surprise. On the other hand, when the initial energy of the lattice is 
larger then a certain threshold, equipartition indeed occurs. 

For a theoretical understanding of the Fermi-Pasta-Ulam experiment, one has 
always tried to link the FPU lattice to a completely integrable system. These are 
dynamical systems possessing a complete set of integrals of motion and therefore 
they display the regular type of behavior that was observed in the FPU experiment. 
With this in mind, two main approaches have been proposed. Firstly, some unex- 
pectedly regular phenomena were observed numerically by Zabuski and Kruskal 23 
in the Korteweg-de Vries (KdV) equation. It was later proved by Gardner et al. 
[Ej that the latter equation is integrable. In fact, Peter Lax realised that KdV is a 
member of a hierarchy of integrable equations that have a Lax-pair representation 
and therefore a complete set of integrals. See []"""] f° r a good overview of these 
results. Using the line of thought of Lax, Flashka 8 j could for instance prove the 
complete integrability of the so-called Toda-lattice, which is a finite-dimensional 
Hamiltonian system very similar to the FPU lattice. On the other hand, one may 
hope to construct an integrable partial differential equation describing the behav- 
ior of the FPU lattice with some accuracy. Thus one assumes the existence of a 
smooth interpolation function u : [0, 1] — ► R, and writes qj = u(j/n). If we now let 
n grow to infinity and choose appropriate ra-dependent scalings of the Hamiltonian 
function of the FPU system, we can -at least formally- obtain a partial differential 
equation for u and hope to find that it is an integrable one. It is important to realise 
though that this procedure is not so well-defined: different ways of approximating 
the derivatives of u may lead to different partial differential equations for u, see 
also [H] in which among others a Boussinesq equation is obtained. Moreover, it is 
a priori not clear if the solutions of the resulting PDE constitute good approxima- 
tions for q, even on finite time intervals. Although several heuristic arguments are 
available that link FPU to for instance KdV, see again [TI], I have not been able 
to find any proof of such a statement. 

The second approach differs greatly from that of the continuum approximations 
and is based on finite dimensional considerations -we shall pursue this approach in 
the rest of this paper. As is well-known 2], periodic and quasi-periodic motion is 
typical in completely integrable finite-dimensional Hamiltonian systems. Unfortu- 
nately, the finite dimensional FPU lattice does not belong to this category. One 
possible explanation of the recurrent behaviour of the lattice is therefore based on 
the famous Kolmogorov-Arnol'd-Moser (KAM) theorem [2], [""]. This theorem ex- 
plains that large measure Cantor sets of quasi-periodic motions can also exist in 
classes of nonintegrable Hamiltonian systems, namely small perturbations of cer- 
tain integrable Hamiltonian systems. The only restrictive requirement is that the 
integrable system that we are perturbing satisfies a certain nondegeneracy condi- 
tion, which requires that each quasi-periodic motion of the integrable system has a 
different frequency. Even though various -again heuristic- arguments advocate this 
approach, and I mention in particular [T""", the big problem is that it is not at all 
a priori clear whether the FPU lattice can really be viewed as a perturbation of 
such a nondegenerate integrable Hamiltonian system. The only obvious integrable 
approximation to the FPU lattice is its linearisation, which is highly degenerate as 
its frequency map is constant. This problem was pointed out again recently in the 
excellent review paper by Ford [""] and the book by Weissert [2*2*] . 
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An interesting attempt to prove the applicability of the KAM theorem was 
made by T. Nishida ^3], who in 1971 published a paper that considers the FPU 
lattice with a finite number of particles, fixed endpoints and symmetric potential 
energy density function (the so-called /3-lattice). Assuming a rather strong nonres- 
onance condition on the frequencies of this lattice, Nishida computes its so-called 
Birkhoff normal form and shows that this normal form constitutes a nondegenerate 
integrable approximation to the original lattice Hamiltonian. Thus he proves the 
applicability of the KAM theorem and the existence of a positive measure set of 
quasi-periodic motions in the nonlinear FPU lattice. But all of this is under the 
assumption of a nonresonance condition, which unfortunately is only satisfied in 
exceptional cases. The actual value of Nishida's computation therefore remains 
unclear. 

This paper is nevertheless devoted to a full proof of what Nishida intended to 
show. The main result can be summarized as follows: 

The Fermi-Pasta- Ulam lattice with fixed endpoints and an arbitrary finite number of 
moving particles possesses a completely integrable finite order Birkhoff normal form, 
which constitutes an integrable appoximation to the original Hamiltonian function. 
The integrals are the linear energies of the Fourier modes. When the potential 
energy density function of the lattice is an even function ((3-lattice), this integrable 
approximation is nondegenerate in the sense of the KAM-theorem. This proves the 
existence of a large-measure set of quasi-periodic motions in the low-energy domain 
of the (3-lattice. 

The key to proving this result lies in the fact that Nishida's nonresonance condition, 
which a priori seems highly necessary for computing the Birkhoff normal form, is 
actually obsolete. As in ^H], which treat the FPU lattice with periodic boun- 
dary conditions, discrete symmetries are the key to proving Nishida's 'conjecture'. 
The results of the present paper can be considered as an extension of ^SI to the 
lattice with fixed endpoints with a considerably simpler proof which again uses 
discrete symmetry together with a simple algebraic trick. 

I want to remark here that the results of this paper do not provide any explicit 
bounds on the domain of validity of the normal form approximation. In particu- 
lar we have at this moment no estimates on the behaviour of this domain when n 
grows to infinity. The principal interest of the result lies in the fact that, at least to 
my knowledge, it is the first complete proof of the very existence of quasi-periodic 
motion in the FPU lattice with fixed endpoints. 

2 The lattice equations of motion 

As was mentioned before, the physical variables of the FPU lattice are the positions 
and conjugate momenta (q,p) = (qi, . . . , q n ]Pi, ■ ■ ■ ,Pn) of the particles in the lattice, 
of which we assume here that there are only finitely many. These positions and 
momenta are elements of the 2n-dimensional cotangent bundle T*R n = R 2n of 
R n . T*R n is a symplectic manifold with the canonical symplectic form dq A dp := 
X^jLi dqj A dpj. Given a Hamiltonian function H : T*R n — * R, the Hamiltonian 
vector field Xh on T*R n is implicitly defined by the relation dq A dp(Xn, •) = dH. 
The integral curves of Xh therefore are the solutions of the system of ordinary 
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differential equations 

, _ dH . dH 

^ dpj ' ^ dqj 

In the case of the FPU lattice the Hamiltonian function is the sum of the kinetic 
energies of all the particles and the interparticle potential energies: 

3 

in which W : M. — > R is traditionally a potential energy density function of the form 

W(x) = \x 2 + |x 3 + ^ (2.2) 

The parameters a and (3 measure the nonlinearities in the forces between the par- 
ticles in the lattice. 

The range of the summation index j in (|2,1|) depends on the exact boundary 
conditions that we impose on the lattice. For lattices with fixed endpoints the mo- 
ving particles are labeled j = 1, . . . , n and the endpoint particles are kept at rest, 
i.e. qo = q n +i = for all time. When we impose periodic boundary conditions, we 
label the particles by elements of the cyclic group, j E Z/ tvz, so that the first and 
the last particle are identified, that is qo = qjy for all time. For conveniency we 
have denoted the number of particles in a periodic lattice by N (i.e. not by n), so 
that its phase space becomes T*~R N . 



3 Discrete symmetry 

The Hamiltonian function (|2.1|) of the periodic FPU lattice (i.e. summation over 
j G "L/ ni) has discrete symmetries of which we shall discuss some dynamical con- 
sequences. Two important symmetries of the periodic FPU lattice are the linear 
mappings R, S : T*R N -> T*!^ defined by 

R :(qi,q 2 , ■ ■ ■ , qN-1, Qn;Pi,P2, ■ ■ ■ ,PN-1,Pn)^ 

(<?2,<73, • • • ,qN,qi\P2,P3, ■ ■ ■ i Pn > Pi ) 
S :(<?i, <?2, • • ■ ,qN-i,qN]Pi,P2, ■ ■ ■ ,Pn-i,Pn)^ 

- (qN-i,qN-2, ■ ■ ■ ,qi,qN;pN-i,PN-2, ■ ■ ■ ,pi,pn) 

It is easily checked that R and S are canonical transformations that leave the pe- 
riodic FPU Hamiltonian (|2.1|) invariant, i.e. R*(dq A dp) = S*(dq A dp) = dq A dp 
and R*H{= HoR) = S*H{= HoS) = H. This implies that R*X H = X R * H = X H 
and S*Xh = X$*h = Xh, that is R and S conjugate the Hamiltonian vector 
field Xh to itself. This in turn implies that R and S commute with the time-t 
flows e tXn of Xh- Canonical diffeomorphisms with this property are called sym- 
metries of H and the group of symmetries of H is denoted Gh- The subgroup 
(R, S) = {Id, R,R 2 ,..., R?*- 1 , S, SR, SR 2 , SR N ~ 1 } C G H is isomorphic to the 
iV-th dihedral group, the symmetry group of the A^-gon, as its elements satisfy 
the multiplication relations R N = S 2 = Id, SR = R~ 1 S. As R and S are linear 
mappings, the elements of (R, S) actually define a representation of D?q in T*~R N 
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by symplectic mappings. 

For every subgroup G C Gh, we define the fixed point set 

FixG = {(q,p) G T*R N \P(q,p) = (q,p) VP G G} (3.1) 

Let (q,p) G Fix G and P G G. Then P(e tXn (q,p)) = e tXn (P(q,p)) = e tXlI (q,p), 
i.e. e tXH (q,p) G FixG. Thus we see that FixG is an invariant manifold for the flow 
oiX H . 

Proposition 3.1 When G is compact and consists of linear symplectic isomor- 
phisms ofT*R n , then FixG is a symplectic manifold with the restriction to FixG 
of dq A dp as symplectic form. This implies that whenever Xjj is tangent to Fix G, 
in particular when H is G-symmetric, 

(Xh)\fixG = X (H\ FixG ) 

Proof: Clearly, FixG = rip g Gker(P — Id) is a subspace of T*R N , and hence a 
submanifold. It remains to be proven that for every (q,p) G FixG, the restriction 
of dq A dp to TL jP )(Fix G) C XL jP ) (T*R ) is nondegenerate. Let us first of all 
identify T^ p) (T*R N ) by T*R N and (Fix G) by FixG and moreover note that 
since G is compact, it contains a unique left-invariant probability measure l dP\ 
called the Haar-measure of G. We can therefore define the operator 

av G : T*R N -» FixG , v t-* [ P{v)dP 

Jg 

The operator av G is a projection of T*R N onto FixG, called 'averaging over G'. 
Now let v G FixG and w G T*R . Then one easily computes that 

(dq Adp)(v,w G (w)) = (dq Adp)(v, \ P(w)dP) = / (dq A dp)(v, P(w))dP 

Jg Jg 

= \ (dq Adp)(P(v),P(w))dP = / (dq A dp)(v,w)dP = (dq A dp)(v,w) 
Jg Jg 

where in the second equality we have used the linearity of dq A dp in its second 

argument, in the third equality the fact that v G FixG, and in the fourth equality 

that every P G G is symplectic. We now observe that when v G Fix G and (dq A 

dp)(v,w) = for every w G Fix G, then (dq A dp)(v,w) = even for every w G 

T*R N . Hence dqAdp, when restricted to FixG = T( 9p )(FixG), is a nondegenerate 

anti-symmetric bilinear form. In other words, Fix G is a symplectic subspace of 

T*M. N . The final statement of this proposition follows trivially from this result. □ 

Let us now look at a periodic FPU lattice with an even number N = 2n + 2 of 
particles. Then an invariant subsystem is formed by the fixed point set of the 
compact group (S) = {Id, S}: 

Fix (S) = {(q,p) G T*R N \qj = -q2n+2~j,Pj = ~P2n+2-j Vj} 

Clearly, in Fix (S), qo = q n +i = Po = Pn+i = 0. Thus we see that Fix (S) is filled 
with solutions (qi (t) , . . . , q N (t) ; p x (t) , . . . , p N (t)) for which the (q% (t), . . . , q n (t) ; p\ (t) , 
. . . ,p n (t)) constitute the general solution curves of the FPU lattice with fixed end- 
points and n moving particles. Hence, the FPU lattice with fixed endpoints and n 
particles is embedded in the periodic lattice with 2n + 2 particles. By Proposition 
13.11 it can be described as a Hamiltonian system on Fix (£), which has the restric- 
tion of dq A dp as symplectic form, and is determined by the Hamiltonian function 
F ix As coordinates on FixS* one could choose (q±, . . . , q n ;pi, ■ ■ ■ ,p n )- 
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4 Quasi-particles 

Of course, the representation of Dn on T*M. N is the sum of irreducible representa- 
tions. It is quite natural to choose coordinates on T*W N that are adapted to these 
irreducible representations. For the periodic lattice, we thus make the following 
real-valued Fourier transformation. For 1 < k < y define: 

IT ^ . 2jk7r IT ^ . 2jkn 

AT ^ 2jfcvr /T ^ 2jfc7r 



and if A" is even: 



A" ' ViV 



AT , ~ V A r 

jez/ivz iez/jvz 

The new coordinates (Q, P) are called quasi-particles or phonons. The transfor- 
mation (g,p) i-> (Q,P),T*R N T*]^ is symplectic and one can express the 
Hamiltonian in terms of Q and P. If we write for 1)2. ljl 

-ff = i?2 + H3 + i?4 

where P2 is a quadratic polynomial in (q,p) and P3 and P4 cubic and quartic 
polynomials in q, then we find that (see [Tl], |17j or |2U| ^ 



- 1 

^2 = ^2^+^) C 4 - 1 ) 



fc=l 

in which for k = 1, . . . , N the numbers tu^ are the well-known normal mode fre- 
quencies of the periodic FPU lattice: 

uj k : = 2sin(^) 

This means that written down in quasi-particles, the equations of motion of the 
harmonic lattice (a = (3 = 0) are simply the equations for A^— 1 uncoupled harmonic 
oscillators and, as ujm = 0, one free particle. In fact, the Hamiltonian system is 
Liouville integrable in this situation. Integrals are for instance the linear energies 

E k := l -{P% + ulQl) 

The FPU model is of course much more interesting when the forces between the 
particles are nonlinear, i.e. when a or f3 is nonzero. The normal modes then interact 
in a complicated manner that is governed by the Hamiltonians H r (r = 3, 4), which 
are of the form 

JV-l 

0:|0|=r fe=l 
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Here the 9 are multi-indices and the eg are real coefficients. Note that for every 

value of q and 3, H is independent of Qn = —m JZnQi- Hence the total momen- 

v N — J J 

turn Pjy = -^7= Pj is a constant of motion and the equations for the remaining 
variables are completely independent of (Qn,Pn)- It is common to set the latter 
coordinates equal to zero, thus remaining with a system on T*M W_1 with coordi- 
nates (Qi, . . . , Qn-i, Pl> ■ ■ ■ j P/v-i)- As uj\, . . . , usn-i > 0, we can conclude by the 
Morse-Lemma or Dirichlet's theorem pQ, that the origin (Q,P) = is a dynami- 
cally stable equilibrium of this reduced system. 

Assume again that N = 2n + 2. From the definition of the quasi-particles and the 
definition of S, we conclude that S acts as follows in Fourier coordinates 

S : (Qi,... ,Qn-i;Pi,...,Pn-i) i-> 

(Ql, • • • j Qn, —Qn+l, ■ ■ ■ , —QN-l',Pl, ■ ■ ■ , Pn, ~ Pn+1, ■ ■ ■ j — P/V-l) 

So that 

Fix(S) = {(Q, P) 6 T*!^ 1 | Q fc = P fc = 0Vn+l</c<iV-l} 

which is a symplectic manifold isomorphic to T*W n . Using the coordinates (Qi, ■ ■ ■ , 
Q n ; Pi, ... , P n ) on Fix(S'), the restriction of the symplectic form Y^,f=i dQjAdPj to 
Fix (S) is simply X)j=i ^Qj ^dPj. By Proposition |23 the Hamiltonian of the fixed 
endpoint lattice thus is simply the restriction of the periodic FPU Hamiltonian 1)4. H 
IO) to Fix(S), that is 

n 1 

#lFix<5> =^2 (Pfe2 + ^ g '' ) + jff3(gi '---' gn ' '---' 0) + F4(Ql '---' Qn '°'---' 0) 
fc=i 

To distinguish we have used the notation Q k := = 2sin( 2 ^ 2 ) (1 < k < n) for 
the linear frequencies of the fixed endpoint lattice. 



5 The Birkhoff normal form 

Nishida's idea was to study the Hamiltonian of the fixed endpoint lattice using 
Birkhoff normalisation, which is a way of constructing a symplectic near-identity 
transformation of the phase-space with the purpose of approximating the original 
Hamiltonian system by a simpler one. The study of this 'Birkhoff normal form' 
can lead to important conclusions about the original system. For r > 2, let J- r 
be the finite-dimensional space of homogeneous r-th degree polynomials in (Q, P) 
on T*!^ 1 and let T := © r > a< F r - With the Poisson bracket {•,•}: T x T -> T 
defined by 

V—— - dF dG 

~ dq k dp k dp k dq k 

T is a so-called graded Lie-algebra, which means that {T r ,Ts\ C -P r +s-2- With 
this definition, we have for each F G J 7 , the 'adjoint' linear operator 

ad F :F^F , G^{F,G} 

We recall the following result, a complete proof of which can be found for instance 
in [U, and [TU]. 
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Theorem 5.1 (Birkhoff normal form theorem) Let H = H 2 + H 3 + . . . £ T 

be a Hamiltonian on T*~R N ~ 1 such that H r £ J- r for each r and 

adn 2 ■ G h-» {H 2 , G} , T r — > J> 

is semi-simple (i.e. complex- diagonalizable) for every r. Then for every finite s > 3 
i/iere is an open neighbourhood £ U C T*M. N ~ 1 and a symplectic diffeomorphism 
$ : U -> T^R^" 1 with the properties that $(0) = 0, Z>$(0) = Id and 

<&*tf = tf 2 + #3 + • • • + + 0(\\(Q, P)\\ s+l ) 

where 

ad H2 (H r ) = 

for every 3 < r < s. T/ie transformed and truncated Hamiltonian H := H 2 + + 
. . . + -ff s is called a Birkhoff normal form of H of order s. 

Idea of proof: For H,F £ T, the curve t i-> (e tXp )*H = H o e* Xj? in J" satisfies 
the linear differential equation and initial condition 

j t (e tXF YH = -ad F ((e tXF )*H) , {e 0XF )*H = H 

This implies that 

(e- XF )*H = e* dF (H) = H + {F,H} + ^{F, {F, H}} + . . . 

The transformation $ is now constructed as the composition of a sequence of time- 
— 1 flows e~ XFr (3 < r < s) of Hamiltonian vector fields Xp r with F r £ J- r . The 
idea is that we first choose F 3 £ !F 3 , such that is transformed into 




When ad# 2 is semi-simple, then T 3 = ker ad# 2 © im ad# 2 and we can decompose 
H 3 = H^ er + H 1 ™ for uniquely determined H^ er £ ker ad# 2 and H 3 m £ im ad# 2 . 
If we now choose F 3 such that ad# 2 (-F3) = H^ m , which obviously is possible, then 
(e- x ?s)*H = H 2 + H 3 + . .. for H 3 = H 3 + {F 3 ,H 2 } = H 3 - a d H2 (F 3 ) = H 3 -H™ = 
H^ er £ ker ad# 2 , i.e. adn 2 {H 3 ) = 0. We continue by choosing F4 £ T 4 such that 
(e- XF i)*{(e- XF 3)*H) = H 2 + H 3 + H 4 + ... for which ad^ 2 (i7 4 ) = 0, etc. After 
s — 2 steps we obtain H with the desired properties. □ 

The normal form H is usually simpler than the original H because it Poisson com- 
mutes with the quadratic Hamiltonian H 2 . This firstly means that H 2 is a constant 
of motion for H and secondly that the flow t \— > e tXlI 2 is a continuous symmetry of 
H. 

Also, H and H are symplectically equivalent modulo a small perturbation of 
order C(||(Q, P)|| s+1 ). Studying H instead of H thus means neglecting this per- 
turbation term. So we make an approximation error, but this error is very small in 
the low energy domain, that is for small \ \{Q, P)\\- With Gronwall's lemma, precise 
error estimates can be made. 

Finally, I would like to mention the ill-known bijective correspondence between 
the relative equilibria of the Birkhoff normal form and the bifurcation equations for 
periodic solutions obtained by Lyapunov-Schmidt reduction, as is explained in [H]. 

For Hamiltonian systems with symmetry, the following elegant and well-known 
result is often useful, see [I] and [TU] : 
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Theorem 5.2 Let H = H2 + H3 + ... G J- and G be a group of linear symplectic 
symmetries of H . Then a normal form H = H2 + H3 + . . . + H s for H can be 
constructed such that also H is G-symmetric. 

This is obvious when one realizes that the 'generating functions' F r of the proof of 
Theorem 15. II can be chosen G-symmetric as well. 

We shall also use the following result on normal forms of symmetric subsystems, 
which trivially follows from Proposition 13.11 and the proof of Theorem 15.21 as the 
transformations e~ Fr induced by symmetric Hamiltonian functions F r leave Fix G 
invariant. 

Corollary 5.3 Let H be a Hamiltonian function with compact symmetry group 
G consisting of linear symplectic mappings. Then the normal form of H\pi x G *s 
simply the restriction of the symmetric normal form H of H to FixG, i.e. 

-Hi Fix G = -Hi Fix G 



This corollary tells us that it is sufficient to compute the normal form of the full 
system to know the normal forms of its symmetric subsystems. In particular, to 
find the normal form of an FPU lattice with fixed endpoints, it suffices to know the 
normal form of the appropriate periodic lattice. Normal forms of periodic lattices 
have been studied elaborately in ^ 



6 Nishida's conjecture 

In his 1971 paper, Nishida proved the following result: 

Theorem 6.1 (Proven by Nishida in |15j) Consider the FPU lattice with fixed 
endpoints, a = 0, (3 7^ and n arbitrary. Assume moreover the fourth order 
nonresonance condition on the Qj. = 2 sin( 2 ^ 2 ) (1 < k < n) requiring that 

n n n 

J ^2(lk—m^)Qk / 0VI,m G {0,1,2,...}" with = 4 and |^— m&| 7^ 

k=l k=l k=l 

Then the quartic Birkhoff normal form H = H2 + H4 of the lattice is a function 
of the action variables := Ek/Vt^ (1 < k < n) only and is therefore integrable. 
Moreover it satisfies the Kolmogorov nondegeneracy condition 

det/^/0 
da k da k > 

This implies that almost all low-energy solutions of the (5 -lattice with fixed endpoints 
are quasi-periodic and move on invariant tori. More precisely, the relative Lebesgue 
measure of all these tori lying inside the small ball {0 < H < e}, goes to 1 as e 
goes to 0. 

As we shall see later, the numbers 

n n 

y^fa - m k )Vt k , for ^ \ l k\ + \ m k\ = 4 
k=l k=l 
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are simply the eigenvalues of ad# 2 on T\. Nishida's requirement that they be 
nonzero except in the trivial case that If. = for all k thus just means that the 
subspace ker ad# 2 6 T\ in which H4 must lie is very low-dimensional. It must 
therefore be remarked here that the integrability of the normal form follows almost 
trivially from Nishida's nonresonance assumption. Nishida's article consists mainly 
of the explicit computation of the normal form H of H under the nonresonance 
assumption in order to check its nondegeneracy. 

But unfortunately, resonances do occur, implying that Nishida's nonresonance 
condition is often violated. We have for instance the relations 

sin(vr/6) + sin(37r/14) - sin(vr/14) - sin(5vr/14) = 

sin(vr/6) + sin(137r/30) - sin(7vr/30) - sin(3vr/10) = 

sin(yr/2) + sin(vr/10) - sin(vr/6) - sin(3vr/10) = 

which lead to a violation of Nishida's nonresonance condition if n + 1 is a multiple 
of 21 or 15. 

Nishida refers to an unpublished result of Izumi proving a much stronger non- 
resonance condition on the i n special cases. The result states that no Z-linear 
relations between the £lk exist if n + 1 is a prime number or a power of 2. I was 
not able to trace back Izumi's proof of this statement, but note that a more general 
result had already been obtained in 1959 by Hemmer |12j . who actually derived an 
expression for the total number of independent Z-linear relations between the Qfc 
(1 < k < n) in terms of Euler's phi-function. It turns out that no Z-linear relations 
exist if and only if n + 1 is a prime number or a power of 2. 

Moreover, as the above examples illustrate, resonance relations between 4 eigen- 
values exist for several n and Nishida's condition is therefore sometimes violated. In 
this paper we will nevertheless prove 'Nishida's conjecture' that his theorem holds 
without having to impose any nonresonance condition. 



7 Near-integrability 



Let us start with a review of some observation in [TH] for the periodic FPU lattice. 
First of all we note that, as the symmetry R is symplectic, 

(R* o ad H2 )(G) = R*{H 2 ,G} = {R*H 2 , R*G} = {H 2 , R*G} = (ad Ha ° R*)(G) 

where we have used that H 2 is i?-symmetric. From this result we read off that 
R* and adn 2 commute as linear operators J- r — ► T r . This means that they can 
be diagonalized simultaneously. In |18j this is done by introducing new canonical 
coordinates (Q,P) 1— * (z,Q as follows. For 1 < k < y, we define: 

z k = ^{PN-k ~ iPk) + ^-(Qk + iQN-k) 
ZN-k = —-^{PN~k — iPk) + ~^-(Qk + iQN-k) 

Ck = 2^^ k ~ ^N-k) — -^{QN-k + iQk) 
CN-k = -^—{Pk ~ iPN-k) + \ {QN-k + iQk) 
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and if N is even: 



z f = 7f (g f-^f ) ' c f = 7! (jP f ~ 2iQ f- 



It is then not hard to compute that 



i?2 = / , iUk(zk(k — ZN-k(N-k) + iiON_ZN_C,N_ 



l<fe<f 



which implies that if 0, G {0, 1,2,.. .} w 1 are multi-indices, then 
in which v is defined as 



1/(6,0) := V «j fc (0 fc - W -fc - ©fc + ®N-k) +wn{9n -On) (7.1) 

L — ' 2 2 2 

l<fc<f 

In other words, ad# 2 is diagonal with respect to the basis of J- r consisting of the 
monomials 2®£ e for which |0| + |0| := X^j=i X + l^il = r an d the corresponding 
eigenvalues are the u(Q,6). In particular we observe that ad# 2 is semi-simple on 
every !F r , so that Theorem 15.11 indeed applies. A Z-linear relation in the frequencies 
oo k is called a resonance. For this reason, the monomials z £ for which z/(0, 0) = 
are called resonant monomials. They are particularly important because they are 
exactly the ones that are not in im ad# 2 and thus, as is clear from the proof of 
Theorem 15 .1\ the ones that cannot be transformed away by Birkhoff normalisation. 
As Q k = uj k (l < k < n), Nishida's nonresonance condition is a consequence of its 
analogon for periodic lattices, that can be formulated as follows: 



and 6 k — 0jv-fc — ©fc + &N-k = for each 1 < k < A 



When lei + 101 = 4 and u(Q, 0) = then 9n-Qn=0 

2 2 

2 • 

Of course, this condition is not valid either. 

On the other hand, one may compute, see ^B], that the operator R* : G i— ► Goi? 
acts as follows on the coordinate function z k ,Q k : 

2iTik 2iTik 2iTik 2-7rifc 

R : z k ^ e n Zk , (k ^ e n ( k: z N _ k h-> e n z N _ k , ( N ^ k i-> e jv , 
zjv i— ► — zjv and CiY. I— ► — Cjv 

2 2 2 2 

And as a result we conclude that, as promised, R* acts diagonally with respect to 
the monomials z®C 9 as well: 



in which /i is defined as: 



N 

M (0,0):= >: j(0fc + ©7V-fc-0fc-0^) + -(©^ -On) mo&N (7.2) 

Z 2 2 



2 



By Theorem 15.21 we now know that the normal form of the periodic FPU Hamilto- 
nian must be a linear combination of monomials z®(, 9 that are both resonant and 
symmetric, i.e. for which v(®, 0) = and /i(@, 0) = mod N. The following 
theorem was proven in |18j . The proof below is considerably simpler though. 
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Theorem 7.1 

i) The set of multi-indices (0,0) £ ^y^ 1 f or which |0| + =3, ji(Q,9) = mod N 
and z^(0, 0) = is empty. 

ii) The set of 'multi-indices (0 , 9) G ^^q -1 /or which |0|+|#| = 4, /i(0,#) = mod N 
and v(Q,9) = is contained in the set defined by the relations 9 k — 9^-k ~ @fc + 
0/v— k = 9 n_ — A[ = 0. 

2 2 

Proof: 

Suppose that |0| + \9\ = 3 and /i(0, #) = mod iV. Then we can conclude from 
looking closely at (|7.2|) and (|7.1j) that there must be integers k,l,m ^ mod N 
with/c + / + m = mod N for which i/(6, 0) = 2» sin(^) + 2i sin(^) + 2i sin(^) = 
2zsin(^) + 2isin(^) - 2isin(^ + ^). Now I learnt the following trick from 
Frits Beukers: write 2isin(^) = x — 1/x and 2isin(^) = y — 1/y for some x,y 
on the complex unit circle. Then v(®,9) = x — 1/x + y — 1/y — xy + 1/xy = 
(1 — x)(l — y)(l — xy)/xy. This is zero only in the trivial cases that x = 1 (k = 
mod N), y = 1 (I = mod N) or xy = 1 (m = mod N). But we already 
knew that k,l,m ^ mod N. The result also follows from the convexity of the 
sine function. 

ii) The proof of the second statement is similar but more remarkable, and based on 
the fact that 2isina + 2isin/3 + 2isin7 — 2isin(a + /3 + 7) = x — 1/x + y — 1/y + 
z — l/z — xyz + 1/xyz = (1 — xy)(l — xz){l — yz)/xyz, which again is zero in trivial 
cases only. □ 

In spite of Theorem I7.1| resonances do exist, as was illustrated by the examples 
in Section A full classification of third and fourth order resonance relations in 
the FPU eigenvalues is given in the Appendix to JH]- Resonance relations lead to 
several nontrivial resonant monomials. But according to Theorem l7.1l we now know 
that these nontrivial resonant monomials are not i2-symmetric and hence cannot 
occur in the normal form of the periodic FPU lattice. As a first result, we immedi- 
ately see now that there are no nonzero elements of Tz that are both resonant and 
i?-symmetric. As a result, H$ = automatically. 

To formulate a result for H±, we need to define the following Hop] '-variables. 
For 1 < k < f , let 

Ofc : = o (P£ + P N~k + ^kQk + ^kQN-k) > °fc := QkPN-k - QN~kPk 

ZLU k 

Ck ■= , — (Pi ~ PN-k + ulQl - uj 2 k Q 2 N _ k ) , d k := — (P k P N _ k + uj 2 k Q k Q N - k ) 

ZlD k LU k 



and if N is even 



aN := (P K +ojnQn) 

2 ZUJ N 2 2 2 



2 



Note that H2 can be expressed as 

H 2 = u k a k 



i<k<f 



We moreover observe that when N = 2n + 2, the identities ow = b k = c k = d k = 

TV 
2 



and a k = E k /Vl k (1 < k < hold on Fix (S), so that our definitions agree with 
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the definition of in Theorem 16.11 The following result was proven in JH] for the 
periodic FPU lattice. The proof consists of a careful analysis of the subspace of 
resonant and (R, S^-symmetric polynomials in T± with the help of Theorem 17. II It 
is not very deep and we do not repeat it here. 

Theorem 7.2 Let H = H 2 + H 3 + H 4 be the periodic FPU Hamiltonian £fl\\JJfy . 
Then there is a unique quartic Birkhoff normal form H = H 2 + H4 of H which is 
(R, S)- symmetric. For this normal form we have H3 = 0, whereas H4 is a linear 
combination of the quartic terms a^a^i , bkbk' (1 < k,k' < y) and if N is even also 
ojvOfc (1 < k < y) and dkdN__ k — CkCN__ k (1 < k < ^). 

Corollary 7.3 (Conjectured by Nishida in [15J) Independent of n,a and (5, 

the quartic Birkhoff normal form H = H 2 + H 4 of the FPU lattice with fixed 
endpoints \4-$ is integrable with integrals E^. 

Proof: By Corollary 15.. 3( the Birkhoff normal form of l|4.3|) is the restriction of the 
Birkhoff normal form of the periodic lattice with N = 2n + 2 particles, to Fix (S). 
But on Fix (S), we have that bk = ct = dk = and = Ek/^lk- So according to 
Theorem E21 Ha is a quadratic function of the Poisson commuting E^. So, clearly, 
is H 2 . □ 

Note that to prove the integrability of the normal form of the fixed endpoint lattice, 
we had to use the symmetry of the periodic lattice in which it is embedded. It must 
also be remarked here that it is very exceptional for a high-dimensional resonant 
Hamiltonian system to have an integrable normal form. 

Let us dwell a little longer on the dynamics of the normal form and consider the 
integral map E : T*R n -> R n that sends (Q, P) i-> (E 1 , . . . , E n ). One checks that 
when Ek > for every k, the derivatives DEk(Q, P) are all linearly independent. 
As the level sets of E are moreover compact, the theorem of Liouville-ArnoFd en- 
sures that for each e = (ei, . . . , e n ) with > for each k, the level set -E _1 ({e}) 
is a smooth n-dimensional torus. 

To compute the flow on these tori, we transform to action-angle coordinates 
(Q,P) 1— ► (a, <f) as follows. Let arg : M 2 \{(0,0)} — > M/271-z be the argument func- 
tion, arg : (r cos r sin <I>) 1— ► $ and define 

fk = arg(P fc , SlkQk) , a k = E k /Vt k = 7^r( P k + ^iQk) > 1 <k<n 

With the formula d arg(x,y) = xd y~y^ ; one can verify that (<p,a) are canonical 

x -\-y 

coordinates: dQ A dP = dip A da. So in these coordinates the equations of motion 

read 

n . _ dH 4 (a) 
ak = , ipk = Uk~\ x 

This simply defines periodic or quasi-periodic motion. Remark: (</>, a) are some- 
times called 'symplectic polar coordinates'. 
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8 Nondegeneracy 



To verify that the normal form H is nondegenerate, we examine the frequency map 
si which assigns to each invariant torus the frequencies of the flow on it: 

il : a I— ► I ill H X ; . . . , Si n + 



aoi oa n 

The nondegeneracy condition of the KAM theorem requires that si be a local dif- 
feomorphism, which is the case if and only if the constant derivative matrix gf fc g^ 4 ; 
is invertible. To check this, we will unfortunately need to compute the Birkhoff 
normal form explicitly, where until now we had been able to avoid this. In the next 
Theorem we shall present the normal form of the FPU Hamiltonian in the case that 
H% = 0, i.e. a = 0. This lattice, that has no cubic terms, is usually referred to as 
the /3-lattice. 

Theorem 8.1 (Conjectured by Nishida in [15J) If a = 0, then a quartic Birk- 
hoff normal form of FPU lattice with fixed endpoints is given by H = Hi + H\, 
where 

\l<k<l<n l<k<n 

'Proof: The computation of the normal form had already been performed by 
Nishida JS] who obtained exactly the above normal form, but under the assump- 
tion that resonant monomials are absent in the lattice Hamiltonian. We now know 
that these monomials can not occur in the Hamiltonian as they are not i?-symmetric 
in the corresponding periodic lattice. Hence Nishida's computation gave the correct 
answer. 

The reader can find similar computations in ^S], [2U] and |5J of the normal 
form of the /3-lattice with periodic boundary conditions. We can therefore obtain 
the result also by substituting = E^/Qk, = = d& = on Fix (S) in the 
normal form of the periodic lattice that was obtained for instance in Theorem 10.1 
in [IB]. □ 

It is now an easy excercise to prove the invertibility of the matrix gf fc g^ 4 ; ■ Its 
nondegeneracy was also checked by Nishida himself by applying elementary row 
and column operations to compute the determinant that turns out to be nonzero. 
Thus we conclude: 

Corollary 8.2 (Conjectured by Nishida in [15J) If a = and (3^0, then the 
integrable quartic Birkhoff normal form H = Hi + of the FPU lattice with fixed 
endpoints \4-$ satisfies the Kolmogorov nondegeneracy condition. Hence almost all 
low-energy solutions of the FPU lattice with fixed endpoints are quasi-periodic and 
move on invariant tori. In fact, the relative measure of all these tori lying inside 
the small ball {0 < H < e}, goes to 1 as e goes to 0. 

Nishida, and we, chose to compute normal form H2 + H4 only for the /3-lattice. 
This computation is already quite long, but it becomes extremely hard when 
It should nevertheless also be possible to write down an expression for the fixed 
endpoints normal form if a 7^ 0. For checking Kolmogorov's condition this will 
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actually be necessary. We know a priori that the resulting normal form will be 
integrable and depends quadratically on the E^ (or a&). It is very likely that for a 
large open set of a and (5 the nondegeneracy condition holds and the KAM theorem 
applies. Without computation this is clear for \a\ <C |/3| (and n fixed) because then 
the coefficients of the normal form can differ only slightly from those in Theorem 

9 Acknowledgement 

The author would like to address many thanks to F. Verhulst, J.J. Duistermaat, 
G. Gallavotti and F. Beukers for many valuable remarks and comments, without 
which this paper would not have been written. 

References 

[1] R. Abraham and J.E. Marsden, Foundations of mechanics, Ben- 
jamin/ Cummings Publishing Co., Reading, Mass., 1978, Second edition. 

[2] V.I. Arnol'd and A. Avez, Ergodic problems of classical mechanics, W. A. 
Benjamin, New York-Amsterdam, 1968. 

[3] H.W. Broer, G.B. Huitema, and M.B. Sevryuk, Quasi-periodic motions in 
families of dynamical systems, Springer- Verlag, Berlin, 1996. 

[4] R.C. Churchill, M. Kummer, and D.L. Rod, On averaging, reduction, and 
symmetry in Hamiltonian systems, J. Differential Equations 49 (1983), no. 3, 
359-414. 

[5] R.H. Cushman and L.M. Bates, Global aspects of classical integrable systems, 
Birkhauser Verlag, Basel, 1997. 

[6] J.J. Duistermaat, Bifurcation of periodic solutions near equilibrium points of 
Hamiltonian systems, Bifurcation theory and applications (Montecatini, 1983), 
LNM 1057, Springer- Verlag, Berlin, 1984, pp. 57-105. 

[7] E. Fermi, J. Pasta, and S. Ulam, Los Alamos report LA-19^0, E. Fermi, Col- 
lected Papers II (1955), Univ. Chicago Press, Chicago, 1965, pp. 977-988. 

[8] H. Flaschka, The Toda Lattice. II. Existence of integrals, Phys. Rev. B 9 
(1974), 1924-1925. 

[9] J. Ford, The Fermi- Pasta- Ulam problem: paradox turns discovery, Phys. Rep. 
213 (1992), no. 5, 271-310. 

[10] G. Gaeta, Poincare normal and renormalized forms, Acta Appl. Math. 70 
(2002), no. 1-3, 113-131, Symmetry and perturbation theory. 

[11] C.S. Gardner, J.M. Greene, M.D. Kruskal, and R.M. Miura, Method for solving 
the Korteweg-de Vries equation, Physics Rev. Lett. 19 (1967), 1095-1097. 

[12] P.C. Hemmer, Dynamic and stochastic types of motion in the linear chain, 
Thesis, University of Trondheim, Trondheim, 1959. 

[13] F.M. Izrailev and B.V. Chirikov, Statistical properties of a nonlinear string, 
Sov. Phys. Dokl. 11 (1966), 30-32. 



15 



[14] E.A. Jackson, Nonlinearity and irreversibility in lattice dynamics, Rocky 
Mountain J. Math. 8 (1978), no. 1-2, 127-196. 

[15] T. Nishida, A note on an existence of conditionally periodic oscillation in a 
one- dimensional anharmonic lattice, Mem. Fac. Engrg. Kyoto Univ. 33 (1971), 
27-34. 

[16] R.S. Palais, The symmetries of solitons, Bull. Amer. Math. Soc. (N.S.) 34 
(1997), no. 4, 339-403. 

[17] P. Poggi and S. Ruffo, Exact solutions in the FPU oscillator chain, Physica D 
103 (1997), no. 1-4, 251-272. 

[18] B. Rink, Symmetry and resonance in periodic FPU chains, Commun. Math. 
Phys. 218 (2001), no. 3, 665-685. 

[19] , Symmetric invariant manifolds in the Fermi- Pasta-Ulam lattice, 

Physica D 175 (2003), 31-42. 

[20] B. Rink and F. Verhulst, Near-integrability of periodic FPU-chains, Physica A 
285 (2000), 467-482. 

[21] J. A. Sanders, On the theory of nonlinear resonance, Thesis, University of 
Utrecht, Utrecht, 1979. 

[22] T.P. Weissert, The genesis of simulation in dynamics, Springer- Verlag, New 
York, 1997. 

[23] N.J. Zabuski and M.D. Kruskal, Interaction of "Solitons" in a collisionless 
plasma and the recurrence of initial states, Phys. Rev. Lett. 15 (1965), 240- 
243. 



16 



